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Abstract 

It is shown that the standard expression for the diffusive flux in 
a semiconductor can be transformed to the flux from the Einstein 
- Smoluchowski equation after making use of the "Einstein identi- 
ties". Later, it is possible to use a transformation introduced in: 
S. Piekarski, "On the modified Fick law and its potential applications" 
(J.Tech.Phys., 44,2, 125-131,2003). In this formalism, it is possible to 
formulate simple models of a semiconductor and an electrolite and to 
obtain the alternative description of a Debye effect (by the Debye ef- 
fect we mean here the screening of the point particle potential). Final 
results are compared with the standard analysis. 



1 Introduction 



According to Peierls ([1], p. 204, Eq. (10.21)), the electric current J in a semi- 
conductor is given by 

j=- a ^t-D— (I) 

dx dx' 

where a is the conductivity, D is the diffusion coefficient, <ft is t ne electric 
potential and n is the concentration of the electrons. The conductivity a and 
the diffusion coefficient D are related by the "Einstein identity" 

D = — - = kTu, (2) 
e n 
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where u = a/ne is independent of the electron density, and is called the 
"mobility" (compare also [5], p. 199). 

From the "Einstein identity" (2) one can determine the conductivity a as 
a function of the temperature and the diffusion coefficient 



a = 



eDn 



After inserting (3) into (1) one arrives at 

eDn d(j) 



j d(f) dn 
dx dx 



kT dx 



ox 



-D 



dn en d(f) 
dx kT dx 



(3) 



(4) 



It can be seen that (4) is a flux from the Einstein - Smoluchowski equation 
(up to minor differencies in the conventions for signs and notations) [2, 3, 8, 9]. 
For the constant temperature field, this flux can be written in the following 
form [2, 3] 



J = —D exp 
Let us show this explicitly 

—D exp 

r 

—D exp 



e0 


d j 


e<p 






— < n exp 
dx { 




\ 



(5) 



e0 




e0 










\ 



e0 




ecj) 


d 




exp 


kT 


dx 



-n — D exp 



e0 
"kT 



d 

n— exp 

dx 



kT 



dx 



Dexp 



e0 




e0 


d 


e0 




exp 




dx 





P d Dne d ^ 
dx kT dx 



d ne d 
d~x n+ lTd-x^ 



J. 



(6) 



It should be stressed again that the temperature T in (5) and (6) is constant. 

On p. 204 of [1] Peierls writes that the concentration of electrons is given 
by the "Boltzmann formula" if the system is in equilibrium in the spatially 
varying electric potential = (x) . This property can be shown immediately 
after putting the flux (5) equal to zero: 



J 



-D exp 



e0 




e0 




_~kf_ 




kT_ 


} 



0. 



(7) 
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Since the diffusion coefficient is assumed to be different from zero, the above 
condition implies that 



d 



e<j) (x) 
kT 



(8) 



and therefore 



n (x) = C exp 



ecf> (x) 
kT 



C>0. 



(9) 



It is worth to discuss also the explicit form of the Einstein - Smoluchowski 
equation for the electrically neutral matter; in [2, 3], the following convention 
has been used: 



d_ 
dt 



p(x,t) = ^D'(T,x)^p(x,t) + 



p(x,t) d 
kT dx 



U(x) 



(10) 



where U (x) is a potential for the external force acting on the diffusing 
molecule, p is the concentration of the diffusing molecules, and D' (T, x) is 
the diffusion coefficient. The equation (10) can be generalized to the equation 
for diffusion of the electrically charged matter after the replacement 



U(x) -^U(x)-q<f)(x,t) 



(11) 



where q is the electric charge of a diffusing molecule and <fi (x, t) denotes the 
electric potential [4]. The flux of diffusing molecules, corresponding to the 
equation 



d 



p(x,t) 



d 



dt dx 
is given by 



i 



d 



D' pOMH 



p(x,t) d 
kT ~dx~ 



[U(x)-q<P(x,t)}\, (12) 



J' = D> (T, x) | |-p (x, t) + P ~^Yx W (X) " q4> (X > t)] }■ (13) 

In turn, the electric flux, corresponding to the flux of molecules (13) , is equal 
to the product of the flux of molecules (13) and the electric charge of a single 
molecule q : 



J' e = qD'(T,x) 



(14) 
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In particular, if the non - electric potential U (x) is constant 



U (x) — Uo — const. (15) 
the electric flux (14) takes the form 

J> = qD' (T, x) | (x, t) + Wo ~ qd> (*, 01 

^( r> ,){^(« > t) + ^M^ [ -^(« >0 ]} = 

?D' (T, x) ^ P (x, t) - q*D> (T, x) t) . (16) 

The electric flux (16) is identical to the electric flux from the Peierls formula 
under the conditions 

f D " ° 1 'x, t) = -oU (x, t) = -q>D> (T, x) 4^1:0 (*, t) (17) 



kT &r v ' ' v ' ' * y ' y kT dx 

and 

-D—n = qD'(T,x)—p(x,t), (18) 

(in (17), the "Einstein identity" (3) has been taken into account). It implies 
that 

n = p 

and 

e = q. 

From the relation 

-D = qD'(T,x) (19) 

one can see that there is a difference in notational conventions for the diffusion 
coefficients. In the Peierls' version, the charge of the carrier is inserted into 
the definition of the diffusion coefficient while in (14) they are explicitly 
separated. It seems useful to compare both fluxes with the balance of electric 
charge, derived from Maxwell equations. In general, the electric flux is a part 
of Maxwell equations 

rotH = j + ^D e (20) 
rotE = -^-B (21) 



4 



divB = 
divD e = p s . 



(22) 
(23) 



After taking a divergence of (20) and the time derivative of (23) one arrives 
at the identities 

d 

= divj + — x divD e (24) 



and 



dt 



—divD e = — p.. 
dt dt Hs 



Comparing of (24) and (25) gives 



d 



df Ps + divj 



0. 



(25) 



(26) 



From the form of (26) one can see that the sign convention applied in (14) 
is opposite to that used in (20) (compare [2, 3]). In the rest of this paper we 
shall use the notational convention used in (14), but without the prime in 
the notation for the diffusion constant. 

In order to discuss the Debye effect in a semiconductor, let us assume 
that the electric charge p s is a sum of the electric charge of the negative 
carriers en(x,t) and the electric charge density of the static "background" 
— ep b . The processes of ionization and recombination are neglected. Since 
we are interested in the Debye effect, we assume that, besides the constant 
background charge, there is a point charge neS (x) at x = Therefore, the 
final expression for p s is 



p s = en (x, t) — ep b + ne5 (x) 
while the explicit form of the diffusive flux in Maxwell equations is 



n(f) 




ri(p 




7kT. 


l{ n 6XP 




\ 



j = —eD (x, T) exp 
Under these assumptions, Maxwell equations become 
rotH = —eD exp 



e0 


d j 


e0 






— < nexp 
ox { 




\ 



d 

rotE = --B, 
dt 



(27) 



(28) 



(29) 



(30) 
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divB = 0, 



(31) 



divD e = en (x, t) — ep b + neS (x) , 
and the balance of the electric charge becomes 



d_ 

dt 



[en (x, t) — ep b + ne5 (x)] + div < —eD exp 



e<p 
kf 



d_ 

dx 



n exp 



kT 



(32) 



d_ 
di 



en (x, t) + div < — eD exp 



e<p 



d_ 

dx 



nexp 



For simplicity, one can assume that 



and 



B = fi p r H, 
B = 0. 



In turn, (30) and (36) imply that 



rotE = 0. 



If 



E = -V0 (x,i) 
then the remaining equations take the form 



= — eD exp 



and 



kT 



e0 




e0 










} 



(33) 

(34) 
(35) 
(36) 

(37) 
(38) 



+ jU £r [-V0] , (39) 



— €oe r V 2 = en (x, t) — ep b + neb (x) . (40) 

In order to describe the process of dynamical screening one has to impose 
the corresponding initial and boundary conditions. Let us assume that the 
support of the diffusion coefficient (or, strictly speaking, the tensor of the 
diffusion coefficients) is restricted to the interior of the sphere with a radius 
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R. The important class of solutions of (39) and (40) consists of stationary 
solutions for which the time derivative is put equal to zero 



e0 




e0 




_~kf_ 




_kT_ 


) 



= — eD exp 

— EQE r V 2 (j) = en (x) — ep b + r/e5 (x) . 
The condition (41) 1 can be solved to give 

e<j) (x 



50 



0, 



n(x) = C exp 



kT 



,C>0 



(41) 
(42) 

(43) 



(compare (7) — (9)). In principle, the value of the positive paremeter C can 
be determined from the initial conditions. 
Now, (43) can be inserted into (42): 



-£q£t-V 2 (x) = eC exp 



(x) 



kT 



epb + Ve5 (x) . 



(44) 



This equation depends explicitly on the three parameters C, p b and rj. The 
fourth parameter is the radius R of the sphere in which the transport pro- 
cesses take place. In the more detailed analysis one should use the dimen- 
sionless variables, taking into accunt the spatial scales of the system (for 
example, the radius R ), in the way similar to that used in the Chapman - 
Enskog expansion, but in our preliminary discussion we simply assume that 



exp 



e0 (x) 
kT 



1 



After inserting (45) into (44) one arrives at 

ecf) (x) 



-£ £ r V (j) (x) = eC 



1 - 



kT 



ecf) (x) 
kT 



- ep b + ne5 (x) . 



(45) 



(46) 



It is well - known that the solution of the modified Poisson equation 

(V 2 -p 2 )(j) = -4ne5 (r) (47) 

is of the form 



<j)(r) 



- exp (— iir) 
r 



(48) 
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Some solutions of (46) can be determined from the solution (48) of the mod- 
ified Poisson equation (47) after the corresponding change of variables. Let 
us describe it in detail. 
Eq.(46) is equivalent to 

V 2 - r^-) (x) = ~[C- Ph ] - (*) • (49) 

Kl cQt r J t()t r EQt r 

Therefore, one can introduce the notation 

2 ^ ^ { \ 

^ kTSQE r 

Now, we would like to change the variables in 

(V 2 - p 2 ) (x) = -— [C - Pb ] - ^-8 (x) , (51) 

SqS t Eq£ t 

transforming it into the modified Poisson equation. To this aim, we insert 

(x) = 0' (x) + a (52) 

into the l.h.s. of (51): 

(V 2 - ^) [0' ( x ) + a] = (V 2 - /i 2 ) 0' {x)-fa = -— [C - Ph \-^-b (x) . 

E§£ r Eq£ t 

(53) 

Eq.(53) is equivalent to 

(V 2 - p 2 ) 0' (x) = p 2 a- — [C- Pb ] - ^-b (x) , (54) 

and the value of a can be determined from the condition that the sum of 
first two terms on the r.h.s. of (54) should vanish. It means that 

p 2 a - — [C - p b ] = (55) 

what after insertion of (50) into (55) gives 

kT 

a = — [C- Pb ]. (56) 
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Therefore, our equation (51) takes the form 



[V 2 - //] 0' (x) 



rje 



£q£ t 



Six). 



(57) 



In order to transform (57) into the "modified Poisson equation" one can make 
the following transformation 



After comparing (58) and (47) one sees that 

e 2 C 



(58) 



4>'{x) = 



rje 



A7re e r x 



cxp 



[kTeoEr 



x 



(59) 



and therefore 



(x) = <j)' (x) + a 



rje 



AtIEqEtX 



cxp 



e 2 C 
kTeo£j 



x 



k-T 



A similar discussion can be made in the case of an electrolyte. 

Again we start with the Maxwell equations but in thie case of an elec- 
trolyte we introduce two different kind of carriers, one charged with a negative 
charge and the second one charged with a positive charge. The density of the 
negative carriers is denoted n (x, t) while the density of the positive carriers 
is denoted p(x,t). We also introduce two different fluxes 



and 



j n = -eD n exp 



j p = eD p exp 



e0 


d j 


e(f> 






— < nexp 
ox I 




} 



e(f> 


d r 


e<j) 


JcT_ 


— <^pexp 

ox y 





(61) 



(62) 



Also in this case we introduce the point charge r]e5 (x) fixed at the point 
x = 0. We also assume that the static background charge is absent in this 
case and the processes of ionization and recombination are absent. Therefore, 
Maxwell equations are in this case 



rotH = —eD n exp 



e0 


d j 


e0 






dx 1 n6XP 
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eD p exp 



ecp 


d j 


ecp 






— <pexp 
ox I 




\ 



(63) 



d 



rotE = --B 
dt 



(64) 

divB = (65) 

divD e = en (x, t) — ep (x, t) + neS (x) . (66) 

Our futher discussion is analogous to the previously discussed case of a semi- 
conductor. We still assume (34), (35), (36) and obtain (37) as a consequence. 
Then, we assume (38) and obtain, as the counterparts of (39) and (40), equa- 
tions 



= —eD n exp 
eD p exp 



e<p 


d r 


ecj) 






— <^ nexp 
ox y 


_kT_ 


\ 



+ 



e0 


d j 


e0 




_kT_ 


— <pexp 
ox I 







d_ 
dt 



and 



(67) 
(68) 



— e e r V 2 (x, t) = en (x, t) — ep (x, t) + neb (x) . 

We discuss the particular solutions of (67) and (68), with the time derivative 
equal to zero and both fluxes equal to zero 



—eD n exp 
eD p exp 



e0 


d j 


e0 






— < nexp 
ox I 


_kT_ 


\ 







e0 


d j 


e<p 




~dx r exp 





0. 



(69) 
(70) 



Similarly as before, we solve (69) and (70) and insert the results into (68), 
arriving at 



— e e r V 2 (x) = eC n exp 



e<p 
kf 



— eC p exp 



kT 



+ r]e5(x) (71) 



this time. Both exponents from (71) are linearized according to the rules 

exp 



ecp 
kf 



1 _ 

kT 



exp 



kT 



1 + 



ecp 
kf 



(72) 
(73) 
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and therefore one arrives at 
— E E r V 2 (x) = eC n 



1 - 



kT 



— eC n 



1 + 



kT 



+ r]e5 (x) . (74) 



This equation is equivalent to 



[C n + C p ] 



E§E r kT 



(x) = 



EqE, 



[C n — C p ] 



rje 



5 (x) . (75) 



Similarly as in the case of a semiconductor (50), one defines the parameter 
ji 2 which is now 

(76) 



" 2 = i^T |C " + C > J 
and after inserting (76) into (75) one arrives at 

[V 2 -/i 2 ]0(x) = --^[C„-C p ] 

CQC r Ef)E r 

After inserting (52) into the l.h.s. of (77) one arrives at 



S(x). 



[V 2 - fi 2 ] <p' (x) - fi 2 a 



[C n — C p 



rje 



Six) 



EqE t £o&r 

The above furmula is an analogy to (53). Now, (78) is equivalent to 



[V 2 - fi 2 } 0' (x) = fi 2 a 



-^-[Cn-Cp]-^-8{x) 

EqE t EqE t 



(77) 



(78) 



(79) 



and the value of a is determined here from the condition that the sum of 
first two terms on the r.h.s. of (79) should vanish. Therefore 

e 



[i 2 a 



EqE, 



[C n — C p 



what after insertion of (76) into (80) gives 

kT [C n — Cp] 



a 



e [C n + C p ] 



(80) 



(81) 



Similarly as in the case of a semiconductor, one can write our equation in 
the form 

[ V2 -"VW = -4*^W (82) 
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and arrive at the solution 



4>{x) 



4lTEoE r X 



exp — 



kT [C n — C p ] \ kT [C n — C p ] 

~[c n + Cp)J x \ ~[c n + c p y 



(83) 



It is worth to compare our approach with that described in Jackson's mono- 
graph [6] on p. 342: 

.."A nonrigorous derivation of the screening effect described above was 
first given by Debye and Huckel in their theory of electrolytes. The basic 
argument is as follows: 

..."A nonrigorous derivation of the screening effect described above was 
first given by Debye and Huckel in their theory of electrolites. The basic 
argument is as follows. Suppose that we have a plasma with a distribution 
of electrons in thermal equilibrium in an electrostatic potential $. Then 
they are distributed according to the Boltzmann factor e~~^ where H is the 
electric Hamiltonian. The spatial density of electrons is therefore 



Now we imagine a test charge Ze placed at the origin in this distribution of 
electrons with its uniform background of positive ions (charge density — en ). 
The resulting potential $ will be determined by Poisson's equation 



n (x) = n exp 



kf 



(84) 



V 2 $ = —AnZeS (x) — Anerio exp 




(85) 




is assumed small, the equation can be linearized: 



V 2 $ - k%$ - -4irZe5 (x) 



(86) 




2 "-0' 



(87) 



...Equation (87) has the spherically symmetric solution: 



$ (r) = Ze 



r 



(88) 



12 



showing that the electrons move in such a way to screen out the Coulomb 
field of a test particle in a distance of the order of k^ 1 ." ... 

It can be seen that the approach described by Jackson contains less pa- 
rameters. The other difference is that in our approach one can state the 
boundary conditions while Jackson considers an infinite system. 

It is also worth to stress that in our approach it is possible to describe 
the relaxation of the system towards the stationary state. 

2 Conclusions 

It would be interesting to solve our equations without linearisations and to 
compare our results with the experimental results (for isothermal conditions). 

In principle, our approach can be applied also to the systems with the 
temperature changing in space and time (compare [2,3,4]). We hope to 
discuss that in future. 

Acknowledgement I am grateful to prof J.Lewandowski and Dr.G.Radzikowska 
for their help. 
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